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We study separability properties in a 5-dimensional set of 
states of quantum systems composed of three subsystems of 
equal but arbitrary finite Hilbert space dimension. These are 
the states, which can be written as linear combinations of per- 
mutation operators, or, equivalently, commute with unitaries 
of the form U ® U ® U . We compute explicitly the following 
subsets: (1) triseparable states, which are convex combina- 
tions of triple tensor products, (2) biseparable states, which 
are separable for a twofold partition of the system, and (3) 
states with positive partial transpose with respect to such a 
partition. 
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I. INTRODUCTION 

One of the difficulties in the theory of entanglement 
is that state spaces are usually fairly high dimensional 
convex sets. Therefore, to explore in detail the potential 
of entangled states one often has to rely on lower dimen- 
sional "laboratories" . An example of this was the role 
played by a one-dimensional family of bipartite states [Q , 
which has come to be known as "Werner states" . In this 
paper we present a similar laboratory, designed for the 
study of entanglement between three subsystems. The 
basic idea is rather similar to Q, and we believe this 
set shares many of the virtues with its bipartite counter- 
part. Firstly, the states have an explicit parametrization 
as linear combinations of permutation operators. This 
is helpful for explicit computations. Secondly, there is 
a "twirl" operation which brings an arbitrary tripartite 
state to this special subset. This proved to be very helpful 
for the discussion of entanglement distillation of bipar- 
tite entanglement: the first useful distillation procedures 
worked by starting with Werner states, applying a suit- 
able distillation operation, and then the twirl projection 
to come back to the simple and well understood subset, 
thus allowing iteration. Geometrically this means that 
the subset we investigate is both a section of the state 
space by a hyperplane and the image of the state space 
under a projection. The basic technique for getting such 
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subsets is averaging over a symmetry group of the en- 
tire state space. Such an averaging projection preserves 
entanglement if it is an average only over local (factoriz- 
ing) unitaries (see [0 for a recent example different from 
ours). 

The third useful property of the states we study is 
that they can be defined for systems of arbitrary finite 
Hilbert space dimension d, which is again important for 
the discussion of distillation. Surprisingly, it even turns 
out that in the parametrization we choose all the sets we 
investigate are also independent of dimension. 

We now describe the entanglement (or separability) 
properties we will chart for these special states. Of 
course, we can split the system into just two subsystems 
and apply the usual separability/entanglement distinc- 
tions. A split 1 1 23 then corresponds to the grouping of 
the Hilbert space H.i®T-L 2 ® H 3 into Hi <g> (H 2 ® H s ). 
We call a density operator p on this Hilbert space 1 1 23 - 
separable (p € B\), or just biseparable if the partition is 
clear from the context, if we can write 



(1) 



(23) 

with A Q > and density operators p a on H.2 ( 3'H3- Fur- 
thermore, as it is a necessary condition for bisepar ability 
(cf. Peres (U), we are going to look at those states having 
a positive partial transpose with regard to such a split. 
Recall that the partial transpose A i— > A Tl of operators 
on Hi <8> H.2 is defined by 



(2) 



where A T on the right hand side is the ordinary trans- 
position of matrices with respect to a fixed basis. 

As a genuinely "tripartite" notion of separability, we 
consider states, called triseparable (or "three-way classi- 
cally correlated" ) , which can be decomposed as 



(i) 

where X a > 0, and the p a are density operators on the 
respective Hilbert spaces. 

The detailed computations leading to the results pre- 
sented here will be published elsewhere, together with 
the discussion of further aspects, such as violations of 
Bell inequalities, and some distillability relations. 
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II. DESCRIPTION OF THE STATES 

We consider only the case, where Hi = TL2 = Ti.3 = TL 
is d-dimensional (2 < d < 00). Then the set of states we 
are going to study, which will be denoted by W, is the set 
of density operators, which commute with all unitaries of 
the form U ®U ®U . The subsets of W we will compute 
will be denoted by B\, for the l|23-biseparable states, V\ 
for the states with positive partial transpose with respect 
to 1 1 23, and by T for the triseparable states. Of course, 
T C B\ cVi C W. We will see that all these inclusions 
are strict, even if we take biseparability with respect to 
all three partitions: T ^ (Bi n B 2 n B 3 ). 

It is a fundamental fact from the representation the- 
ory of classical groups Q that these states are precisely 
those which can be written as a linear combination of 
permutation operators 

peW ^ p = J2 MttKt (4) 

7T 

with coefficients € C and the unitary permutation 
operators V v defined by 

implementing the permutation symmetry of the three 
sites. For density operators hermiticity and normaliza- 
tion reduce the 3! = 6 complex parameters to five real 
ones (an explicit choice will be made below). 
As usual, the integral 

Pp = j dU (U ®U ®U)p{U ®U ®U)*, (5) 

with respect to the normalized invariant measure "dU" 
of the unitary group Ud defines a twirl operation with 
the property p 6 W Pp = p. 

Up to here all statements generalize easily to an arbi- 
trary number of factors, and some are even valid for an 
arbitrary averaging operation with respect to a compact 
symmetry group. However, to carry the analysis further 
one needs at least a precise description of the range of the 
coefficients p, n in equation §Q) , such that the sum indeed 
represents a density operator. In general this is difficult, 
although sometimes (2) one even gets a simplex. In the 
case we study in this paper, the positivity of (^) is best 
seen by using the following basis, rather than the set of 
permutation operators themselves. 

R+= g(l + V( 12 ) + V( 23 ) + V(31) + V (123) + V( 3 21)) (6a) 
R- = g (1 - V (12 ) - V( 23 ) - V(31) + V (123 ) + V (3 21)) (6b) 

i? = I(Vl-F (123) -V (321) ) (6c) 
Ri = l (2 V (23) - V( 31 ) - V (12) ) (6d) 



R2= V3 ( V{12) " V{31) ) (6e) 
Rs = -J= (V (123) - V {321) ) , (6f) 

where we have used cycle notation to represent permu- 
tations. Then , i?_ and Rq are orthogonal projec- 
tions commuting with all V^, and adding up to one. The 
Ri for i = 1,2,3 fulfill the Pauli commutation relations 
[Ri,R±]_ = 0, i2f = R , and i?ii? 2 = iR 3 and cyclic. 
Now every operator p in the linear span of the permuta- 
tions can be decomposed into the orthogonal parts R+p, 
R-P, and RqP, and positivity of p is equivalent to the 
positivity of all three operators. This leads to the follow- 
ing Criterion: 

Criterion 1 For any operator p onTL^TL^TL, define the 
six parameters rk(p) = tr(pRk), for k G {+, — , 0, 1, 2, 3}. 
Then r^iPp) — rk{p). Moreover, each pgWis uniquely 
characterized by the tuple (r + , r_, Tq, r\, r 2 , r 3 ) € R 6 , and 
such a tuple belongs to a density matrix p 6 W if and only 
if 

r+,r_,r o >0, r + + r_ + r = 1 

and r\ + r\ + r\ < rg. (7) 

Taking ro = 1 — r+ — r_ to be redundant, we get a simple 
representation of W as a convex set in 5 dimensions. 

Note that in this parametrization the set W does not 
depend on the dimension d with one exception: for d = 2 
the anti-symmetric projection is simply zero, so for 
qubits we get the additional constraint r_ = 0. Although 
in the case of three Qubits the dimensionality of the set 
is the same as in || the two sets differ. In fact the set 
presented in ]2j can be obtained by averaging over a group 
of local unitaries of order 24. 

With this parametrization we can describe the 5 di- 
mensional set W by the points lying in the (r + , r_ ) plane 
together with a point (r±,r2,rs) in the corresponding 
Bloch sphere of radius ro = 1 — r + — r_. We note that 
a state p G W is invariant under cyclic permutations iff 
n = r 2 = 0, invariant under the interchange 2 «-> 3 iff 
r 2 = ^3 = 0, and invariant under all permutations iff 
n = r 2 = r 3 — 0. The latter set will be denoted by W°. 

III. RESULTS 

The basic results are summarized in Figure 1. On the 
one hand, each point in this triangle corresponds to a 
density matrix p € W with n = r 2 = r% = 0, i.e., a per- 
mutation invariant state. On the other hand, each such 
point stands for the collection of states with the specified 
(r + ,r_), and arbitrary (ri,r 2 ,r3), all of which are pro- 
jected to the state with vanishing (ri, r 2 , r 3 ) upon permu- 
tation averaging. For qubits, one always has r_ = 0, so 
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only the abscissa EB remains. Otherwise all statements 
are valid in any dimension d. 
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FIG. 1. The permutation invariant subset W° C W. Black: 
T n W°; black or dark grey: B\ n W°; light grey: projection 
of B\ onto W° by permutation averaging. Marked points: see 
text. 



A. Triseparable States 

The black triangle ABC in Figure 1 is the set of trisep- 
arable states. Its vertices are obtained by taking suit- 
able pure product states I^X'J'I with 'F = (f>i ® </> 2 ® 03, 
(II 0i II = 1), and applying the averaging projection P. 
Then one gets the points A,B,C, if for i ^ j 



(4>i\4>i) = o 
(4>i\4>i) = i 

(4>i\4>j) = cos(27r/3) 



(8) 
(9) 
(10) 



Note that the "Mercedes Star" configuration for C re- 
quires only two dimensions, whereas A requires d > 3. 
The technique for getting all of T is similar: an arbitrary 
pure product state is projected into W, and T is com- 
puted as the convex hull of these points. This is more 
difficult than it sounds, and the details will be presented 
elsewhere. The result can be summarized as follows. 

Criterion 2 A state p € W is triseparable if and only if 
the following inequalities are satisfied: 



(a) < r_ < i 

(b) \(l - 2r_) < r+ < (1 - 5r_ 



(c) (3rf + [1 



3r+] 2 ) • (1 - 6r_) < 



All extreme points other than A = (i, i, 0, 0, 0) are in 
the plane r_ = 0, i.e., they can also be realized by three 
qubits. For fixed (r + , r_), the shape of T, as embedded in 
the Bloch sphere parametrized by (r\,r2,rz) is shown in 
Figure 2 (innermost convex set). Its threefold symmetry 
is the residue of the permutations of the three sites. 



B. Biseparable States 

We fix the partition 1 1 23 of the system. Since the pro- 
jection of permutation averaging does not preserve bisep- 
arability with respect to this partition, we now have to 
distinguish two sets in Figure 1: those which are bisep- 
arable states with n = r 2 = r 3 = 0, hence permutation 
invariant (represented as dark grey or black), and those 
which are the images of some biseparable state under 
permutation averaging (represented as light grey). So a 
light grey point in Figure 1 has the property that for some 
suitable (ri,r 2 ,r 3 ) one gets a biseparable state. Special 
points with this property are represented by white circles, 
as opposed to filled circles, which lie in the plane. 

One interesting point in Figure 1 is G = (g,0). It is 
biseparable and also permutation invariant. In particu- 
lar, it is biseparable for any partition of the system. But 
it is not triseparable and, in fact, the only extreme point 
of the permutation invariant biseparable set, which is not 
triseparable. It can be obtained by applying P to the 
pure state with vector * = (|112> - j 121) - V3|122))/V5. 

The basic technique for computing B\ is the same as in 
the triseparable case: one takes pure states with vectors 
of the form <F = (pi® (f>23 , and computes the convex hull 
of their images under P. Some special extreme points of 
£>i are given below. They have the additional property 
of being invariant under the exchange of systems 2 and 
3, which is equivalent to r 2 = r% = 0. The following 
table lists the tuples (r + ,r_,ri), and the vector in a 
suitable basis. 



B 
D 

E 
F 



(1,0,0) 

V3' u ' 3> 

(0,0,-1) 

(0,1,-1) 



1111) (11) 

|122) (12) 

(|112)-|121))/V2 (13) 

(|123) - |132))/V2 . (14) 



In addition to these four points there is a sphere 
of extreme points extending also into the r 2 and r 3 - 
dimensions, which is tangent to the line connecting D 
and E. The inequalities describing the biseparable set 
B\ are given in the following 



(ri + r+ - r_) • ((n - 2r+ + 2r_) 2 - 3r$) . 
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Criterion 3 A state p S W is biseparable with respect 
to the partition 1 1 23 if and only if < r_ < |, and one 
of the following conditions holds: 



(a) (3r_ - 1) < (1 + r x - r_ - 2r + ) < ane 



3r 2 



3rd 



(i + 2n 



<(2 



4r 



2r + ) 2 . 



(b) < (1 + n - r_ - 2r + ) < (1 - 3r_) and 

3r| + 3r£ + (1 - 3r_ - 3r + ) 2 < (n + 2r_ - 2r + ) 2 . 



For a typical point in Figure 1, the subset B\ in the 
Bloch sphere is depicted in Figure 2. Note that the 
boundary is composed of two quadratic surfaces, corre- 
sponding to the two alternatives in the above criterion. 



T3 




FIG. 2. The Bloch sphere over the point r+ = 0.27 and 
r_ = 0.1 in W°. Innermost solid: triseparable states T. Be- 
tween T and outer sphere: the boundary of the biseparable 
states 0i. 



invariant subset coincide, so no difference can be seen in 
Figure 1. 

Despite this similarity the technique for computing V\ 
is completely different. The key is the observation that 
partial transposition in the first factor maps operators 
commuting with all unitaries (U ®U ®U) to operators 
commuting with all unitaries of the form (U (8 U ® U). 
Obviously, the latter set is again an algebra, which even 
happens to be isomorphic to the algebra (U <8> U <B> U)- 
invariant operators: two one-dimensional summands plus 
the 2 x 2-matrices. Hence positivity of partial transposes 
can be decided along the same lines as in Criterion 1. 

Criterion 4 Let p G W be a density operator with expec- 
tations rfc = tv(pRk), k = +, — , 1,2,3. Then the partial 
transpose of p with respect to the first tensor factor is 
positive, i.e. pePi, if and only if 



si 



S2 



- n - 

- n - 

and 



5r_ 
r_ 4 



5r_i 



> 

> 



r\ < S1S2/3 



(15a) 
(15b) 
(15c) 



Note that condition (15c) is exactly the same as the 
quadratic inequality in Criterion 3, Part(b). Therefore, 
we need not even provide a new plot for the set V\: In 
Figure 2, this set can be obtained simply by extending 
the quadratic surface, which wraps around T, all the way 
to the boundary of the Bloch sphere. In other words, the 
difference between B\ and V\ is only that states in B\ 
have to satisfy an additional quadratic inequality, which 
is represented in Figure 2 by the surface tangent to the 
Bloch sphere. 
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C. States with Positive Partial Transpose 

Holding the partition 1 1 23 fixed we can compare the 
set of states with positive partial transpose with respect 
to the first subsystem (Pi) to B\ . If Peres' Criterion were 
valid in this case, we would have equality in V\ D B\. It 
turns out that the inclusion is strict, but in several re- 
spects the criterion is amazingly good. To begin with, 
the intersections of both sets with two important hyper- 
planes coincide: namely (1) the r_ = plane (in partic- 
ular, for three qubits), and (2) the r2 = r% — plane, 
i.e. for states which are invariant under the 2 «-> 3 inter- 
change. In particular, the projections to the permutation 
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